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Abstract 

We propose pure quantum parton model of soft interactions. The model is based on 
soft partons which interact eikonally and locally in a transverse plane. We directly con- 
struct S-matrix operator and partons Fock space using developed method of discretization 
of a transverse plane. Counting inelastic states we explicitly derive AGK cancellations 
and cutting rules from relation between inclusive and exclusive parton distributions. We 
also discuss black disk behavior from the point of view of our model. 

1 Introduction 

AGK cutting rules [3] are very important part of high energy phenomenology. Usually it 
formulated using reggeon diagrams. The main problem is the lack of a theory which can 
say how exactly study multipomeron vertexes. Most conventional and naive approaches 
are eikonal and quasi-eikonal models. But there is no theoretical motivation to use this 
models on fundamental level. More predictive and intuitive model can be obtained using 
partons. We can associate multipomeron vertex with existence of multiparton states in 
a hardron wave function. Key task here is to find and define parton properties. Then to 
build the theory using parton as elementary building block. 

There are numerous works about regimes with very high energies and high partons 
density. Formation of a black disk causes especial interest. It has been noticed [3] that 
a t-channel cut of arbitrary reggeon diagram gives only pure n-reggeon exchange and we 
can absorb upper and lower diagram parts to the projectile and target respectively. This 
picture can be naturally reformulated in terms of partons. Hardron is compound object 
builded from partons. If we boost hardron to high Y then in its wave function there will 
be new partons which is strongly correlated with primary partons. During the collision, 
partons evolution is frozen and we deal only with given wave function represented as 
certain vector in Hilbert space. The main observation is that the evolution and scattering 
is distinct and independent parts of the theory. This this similar to the high energy QCD 
where scattering is Wilson lines and evolution is a Hamiltonian-driven task. Obviously, 
the theory of elemental soft partons interactions should be developed. In this paper we 
propose the approach to such model and explicitly derive AGK cutting rules. 
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2 The model 



Our model based on the old Good and Walker idea [T] . Consider interaction of the hardron 
projectile and hardron target at high energy. The projectile, eigenstate of energy, being 
composite, can be described as superposition of states which are eigenstates of quasielastic 
scattering. After scattering, projectile state is not simply proportional to the initial 
projectile. It is superposition of different observable particles. This gives a possibility of 
inelastic diffraction process. In other words, in scattering process coherence of projectile 
lost and disbalance of wave function observed as possibility of decay of the projectile. 

Second ingredient of our model is a partons. In soft hardrons scattering phenomenol- 
ogy a colorless parton models are widely applied. For example, the pomeron is qualitative 
described by ladder-like diagram. It is natural to conjecture that the eigenstates of scat- 
tering which leads to diffraction and partons are closely related. Our hypothesis - that 
they are the same. More precisely, each parton in a projectile wave function is eigenstates 
of quasielastic scattering. Parton can scatter elastically or pure absorptive forming some- 
thing like multiperipheral state where produced particles populate midrapidity region. 
The Fock space of parton states is equivalent to the space of all eigenstates of scattering. 

Third ingredient of our model is a evolution. Instead of studying usual limit Y ^ oo 
we study relatively small y ~ 1 where loop corrections can be neglected. If we increase Y 
then the projectile wave function will change. This happens due to parton emission into 
new opened region in phase space. So instead of increasing Y we can simple change 
projectile wave function and replace old projectile by this new state. The situation is 
very similar to the studying of the high energy scattering in the perturbative QCD where 
JIMWLK evolution equation can be derived by boosting projectile [5]. 

Fourth ingredient of our model is the eikonal approximation for multiparton scat- 
tering. This clause is also QCD-motivated. In QCD amplitude of scattering at mod- 
erate Y can be computed by evaluating Wilson lines of projectile and target partons: 
exp{ig J A^{x)dx'^), where is the classical gauge field generated by all partons [6]. It 
is important that the energetic partons moves along straight lines. From general quan- 
tum field theory framework this picture can be explained as quasiclassical approximation. 
Scattering amplitude can be written as functional integral / e^^^^'^DLp, where (/^(x) - some 
fields related with relevant degrees of freedom, ip must obey appropriate border con- 
ditions which correspond to initial and final energetic partons of projectile and target. 
Quasiclassical solution ipo{x) is a field configuration with extremal action which obey all 
necessary boundary conditions. Scattering amplitude can be computed as e*'^['^°l. 

Fifth ingredient of our model is a transverse locality. It means that classical field 
of each soft parton fast decrease with transverse distance, in opposite to QCD where 
fields decrease polynomially. If the hardron is dilute then the fields of parton in not 
overlap with each other. Geometrically, this means that parton size much smaller than 
average distances between partons and partons can be treated as point-like particles in 
a transverse plane. Hardron transverse extension totally associated with wave function 
of constituent partons, not with parton fields. This is similar to hydrogen atom where 
point-like electron have broad wave function. It should be stressed that high parton 
density is not automatically equivalent to high parton number where our model is not 
applicable. For example, if parton localized in small transverse area then locally density 
is high but there is only one parton. 

Projectile state \^ p) can be written as vector in the Fock space generated by action of 
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parton creation operator ai{x) on the vacuum, where x - transverse position. Respectively 
target state I^'t) generated by action of a^ix). 

Main building block in our model is elementary two parton S'-matrix. Consider colli- 
sion of two partons with positions overlapped in transverse plane. Initial state is 

\^,)=a\al\0) (1) 

Due to all five our assumptions we write 

S\^,)=f,\^,)+,^c^\0) (2) 

where ^,1^ is some constants, and is operator of creation of inelastic state which we will 
call as c-modes. We assume that the created inelastic state is also transverse localized 
and we can create its with local operators c^{x). We can motivate transverse localization 
of c-modes by following reasons. From quasiclassical approximation point of view the c- 
mode is just another quasiclassical solution type of the field equations with special border 
conditions. There are no reasons to expect a discontinuous transverse jump of a classical 
solution. Another argument comes from color flux tube phenomenology. Soft partons 
receives a color charge in collision and color fiux tube between patrons is formed. Such 
tube is also should be transverse localized. 

Constants fi and must obey unitary relation 

+ = 1 (3) 

The next step is to find secondary quantized version of equation ([2]) where arbitrary 
number of partons can present in the initial state. 



3 Discretization of a transverse plane 

Partons Fock space can be constructed by creation operators obeying 

[a{x),a\y)]=5{x-y) (4) 

Although we use boson statistics for partons, our subsequent calculations is not depend 
on parton statistic type. Each parton can have arbitrary transverse position which is 
continuous variable. So this space is not separable. If we deal only with finite partons 
number then the probability to have two partons in some point has zero measure. So the 
occupation numbers can have value only and 1. Note that this is similar to studying 
fermions. 

In order to do more rigorous and simple math we approximate original Fock space by 
Fock space constructed on discrete transverse plane where partons can sits only on lattice 
nodes. Then we can tend lattice spacing to zero. Original parton position we replace 
by nearest node position. In our model if we change lattice size then the the number of 
partons is not change. But in the QCD this is never true due existence of RG group flow 
and DGLAP evolution. Moreover, QCD parton Fock space can be correctly formulated 
only on lattice because if we tend lattice spacing to zero then the parton density grow 
to inflnity. Lattice spacing play role of momentum cutoff. Changing cutoff we go from 
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a coarse-grained description to a fine-grained description. We introduce discrete parton 
creation operators obeying 



[ai,a\]=5ij (5) 



where i,j - indexes of lattice nodes. For arbitrary operator O we denote as : O : normal 
ordered operator, as Oj - its discrete version, as 0{x) - its continuous versions. General 
correspondence rule between continuous and discrete versions of operators is 

Qi = a{xi)V As (6) 

Ni = Nixi)As (7) 

where As - area associated with lattice node and Ni = a|aj operator of parton number. It 
is should be stressed that the ([6]) has meaning only as approximation of matrix elements 
of operators by matrix elements of its discrete version. It is straightforward to show that 

^ iV, = y N{x)dx (8) 

i 

In our model if the lattice spacing is sufficiently small then in each node can be no more 
than one parton. In the continuum limit the probability to have at least two patrons 
with equal coordinates is neglible small. If maximum occupation number is one then we 
have 

:Nf:=0 (9) 

This is very important property for our subsequent calculations. As example we construct 
projector on the Fock vacuum. 

Po = 11(1 - iV,) (10) 

i 

This operator give zero if at least one parton exist. Next we want to find continuous 
version of this operator. With help of Q we can write 

1-N,=: e-^» : (11) 

Po=:e-S^V=:e-/^(")'^" : (12) 
Now we can construct secondary quantized S'-matrix. Formula ([2]) can be rewritten as 

S'i = 1 - (1 - ^)iVi,iA^2,j + i^4"i,'«2,i (13) 

Equation (I13p is very natural and valid only if the lattice spacing greater than parton 
size. Continuous version of full S-matrix can be derived applying derivation method of 
(fT2l) to (pj) 

^ g— (1— /j)As J Ni{x)N2{x)dx+i/^/As J ai{x)a2{x)c^x)dx (14) 

i 

When tending As to zero to get continuous version we have a problem that As can be 
smaller than the parton size sq. Recall that for each As we only have approximation for 
matrix elements and smaller As gives better approximation. So we take As = sq as best 
available approximation. If we go to As < sq then locality in (llSh will lost. When we 
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calculate amplitude of realistic process with given momentum transfer we effectively av- 
erage parton densities over distances of scale of inverse momentum. So it is not necessary 
to have very fine-grained description. Finally, we have 

Q _ g— A J N{x)dx+\ J a{x)c^{x)dx (15) 

where ^ = (1 — /x)sO) ^ = ^^/so and a{x) = ai{x)a2{x), N{x) = Ni{x)N2{x) 

4 Cutting rules 

Now our task is to explore how many c-modes exist in the final state 

I*/) = S\^i) (16) 
We define elastic scattering amplitude as 

M = S\^i) (17) 

m(") = (18) 
where we defined inclusive partons spectrum in the initial state 

/('=)(xi ...Xk) = {^i\N{xi) . . . N{xk)\^i) = /f ^(xi . . . xk)ff\xi ...Xk) (19) 

/(n) = 1. I fin) (^^ . . . ^^^dxi ...dXn (20) 

n\ J 

and where f\ and is projectile and target parton inclusive distributions respectively. 
Note that the factor (—1)"+^ lead to so called shadowing corrections. Definition of the 
total cross section is 

<Jtot = Y.\^n\l-S\^i)\^ (21) 

n 

Using unitary S^S = 1 and complete basis |n) we can obtain optical theorem 

atot = 2ReM (22) 
Inclusive spectrum of the c-modes in the final state is 

/f )(xi ...Xk) = (*/|iVc(xi) . . . N,ixk)\^f) (23) 
where Nc{x) = c^{x)c{x). We define generating functional for inclusive spectrum 

Fc[u] = ^— / /i^^Xl . . . Xn)u{xi) . . . u{Xn)dxi ...dXn (24) 
n=0 ^ 

This functional can be equivalently rewritten as 

FM = {-^f\: e/"(-)^^W'^- : 1^-^) = (^-^iH^ {I + UiN,,i) Si\^i) (25) 
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Using (jl3p and relation (jlip we obtain 

Fc[u] = {^^\ n (l + W\'n,Nl,^N2,^) |^.) = (^.| : g/ l^l'"^-)^^-)'^- : |vl.,) (26) 
i 

and therefore 

/i'=)(xi...x,) = |A|2V('^)(xi...x,) (27) 



Equations (j27p is the direct analog of the well-known AGK cancellations. From (jlSp 
we see that the Z*^*^^ gives the fe-fold parton elastic scattering. Equation ([27|) says that 
A;-fold inclusive inelastic spectrum is solely related to /c-fold contribution to the elastic 
amplitude. 

Now we want to find exclusive distribution of c-modes. In the our model exclusive 
spectrum gives multiplicity distribution in the final state. It can be explicitly constructed 
as 

9n...i„ = (^/l^^c,n . . . iVc,„ n (1 - ^c,)|^/) (28) 

Operator 1 — Nc^i plays role of the vacuum projector at given point x. So we select only 
states with precisely n c-modes placed in nodes ii . . Using definition of inclusive 
distribution (p5]) . equation ([28|) in its continuous version can be expressed as 



5„(X1 ...Xn)- 



6u{xi) . . . 6u{xr, 



(29) 



u=-l 



Hence the generating functional Gc[u] for exclusive distribution can be calculated using 
inclusive distribution 

G,[u] = F,[u - I] (30) 



If we define 



then the normalisation is 



9n = ^ I g^^'Hxi . . . Xn)dxi . . . dXn (31) 



n! 



^g„ = Ge[l] = l (32) 

n=0 

Probability to have n c-modes in the final state equals to g'^. From (j29p we can directly 
calculate integrated exclusive distribution for c-modes 

oo 

9'n = T. 9^n^ (33) 
k=n 

(k) (k) 

where we defined gii as the contribution to g'f^ from fc ■ Now it is useful to relate 

(n) 

probabilities gn to cross sections of subprocesses. Let cj^ is the contribution to the cross 
section of create k number of c-modes from n order. If A; > then we obviously have 

= (35) 
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and for n = we use the following identity 



k 



(T, 



^ a^,^> + + . . . = .^'^^ - .^'^^ + E ,1'^) (36) 

n=0 

From psp it is easy to check that the sum in last expression is zero and we obtain 

-f^=5r+-S=^^r + 2i?eM(^) (37) 
where we used optical theorem ([22]). Substituting fc'^^ from ([27|) into ([55]) we finally have 

= (-l)"+iC,"Z)W 

where 

= (-l)'=+i|Ap'=/('^) (39) 

Equations (|38|) is well known AGK cutting rules (integrated) |3|. We can easy check the 
following property 

E = -loi (40) 

n=0 

From equations (|38p we can derive Pumplin bound for cross section of the diffractive 
dissociation [2j . By definition 

oo 

(^o = ^cri''^ =crel + (^diff (41) 
k=l 

from ([38]) we have 

oo 

Pumplin bound valid only if |A| ^1 and ImA <^ ReA, where A has been defined in (|15|) 
In this limit unitary relation have form 

2ReA = |A|2 (43) 
|^|2fe ^ 2''-^2ReA^ = ]^Re(2Af (44) 



Substituting this into (j42|) . we obtain 

^0 = ^tot - ^2i?e(M/,|l - e-/2^^(-)'^-|M/,) < \jtot (45) 

5 Black disk states 

Note that the eikonal model is equivalent to factorization of inclusive partons distributions 

/(^)(xi...xfc) = /«(xi).../«(xfc) (46) 
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From exclusive distribution point of view this relation can be viewed as possible existence 
of arbitrary number of partons in the wave function. And it is well known that in the 
final state multiplicity have Poisson-like distribution. This situation usually signals about 
some sort of chaos and can arise if statistically many degree of freedom involved in process. 
We only can hope that this scenario will arise in limit of very high Y. But in case of a 
dilute hardron there is no reasons to expect such factorization. Moreover, in case of finite 
number of partons there are no multipomeron exchanges of order higher than partons 
number. 

Consider black disk limit. At first we must clarify meaning of the term "black disk". 
Usually it defined as target with vanishing elastic S-matrix for arbitrary projectile. This is 
rather unpleasant requirement because it is difficult to require exactly = by physical 
reasons. Instead we use more weak version and define black disk as target hardron 
satisfying to relations (j46p with some accuracy. This definition can be used in dilute 
regime. If projectile is the black disk too then from (|14p elastic S-matrix has form 

5M = e-^/^i(^')^2(^)'^^ (47) 

More precisely speaking, our definition of black disk should be called as gray disk. True 
black disk arises at high partons densities. 

Now we show some statistical arguments for inevitability of formation of black disk 
states. At first we consider N partons uniformly distributed in area S in the transverse 
plane. This state can be correctly described by density matrix p. If state have exactly 
partons in the area S then 

p\xi . . . Xn) = ^\xi . . . Xn) (48) 

and it is equal to zero in other cases. Normalization constant is 1/Z = Nl/S^ . Then we 
have 

/(^)(xi) = (iV(xi)) = Sp[N{x,)p] = ^ (49) 
/^^^(-i>-2) = fj(l-^) (50) 

f''Hxi,x„xs) = ^{l + ^-^) (51) 

and so on. We see that if A^ is large than lowest inclusive distributions obey factorizations 
(|46p . So considered state is good approximation to a black disk. Such states can describe 
nucleus wave function in term of partons, at least locally. 

Next argument comes from entropy consideration. Consider large hardron state |^'). 
Let take small area Si within hardron and consider subsystem in Si as quantum system. 
Hilbert space can be factorized as tensor product H' Hi, where Hi - Hilbert space for 
subsystem Si and H' for subsystem S' = S — Si. If we want to calculate averages of 
operators acted in Hi then it is useful to introduce density matrix pi for subsystem 5i. 

= (52) 

where trace taken over H' space. Having non trivial density matrix we can evaluate 
entropy 

Ei = -Spi{plnp) (53) 



8 



where trace taken over Hi space. Now we use general property of chaotic systems that 
entropy tends to be additive. We assume that the hardron is locally chaotic. Equivalently, 
El proportional to Si. Such reasoning similar to derivation of Gibbs distribution from 
microcanonical ensemble. If the entropy of two small subsystems and 5*2 is additive 
then the density matrix of the subsystem 51+2 is a the product pi^2 = Pi ® P2- In 
calculations of f^'^\xi,X2) we can take two subsystems Si and 52 where xi G Si and 

X2 G 52. 

/(2)(a;i,x2) = {N{xi)N{x2)) = Spi{N {xi) pi)Sp2{N {X2) P2) = {N {xi)) {N {X2)) (54) 

So we have factorization (I46p . not for whole hardron but for small areas within it. This 
can be applicable for collision of small dilute projectile with large target. 

6 Discussion 

In this paper we have derived AGK cutting rules in the framework of constructed parton 
model. Remarkable that the result comes from pure algebraic relation between inclusive 
and exclusive distributions for inelastic states. There are some unclarified questions here. 

First question is how inelastic c-modes converts into realistic hardrons. At first sight 
in spirit of color flux tube and multiperipheral states we can view this state as state where 
partons populates whole available rapidity interval but only in same transverse position. 
Then using realistic hardrons wave functions we must project final state on observable 
particles. So partons from different c-modes can compose into one hardron. 

Second question is about nature of our partons. How relate partons to constituent 
quarks or to string modes in dual resonance model? It is not clear. But from many 
examples in field theories we see that excited field modes near some vacuum have Fock 
space form. Crucial assumption is only transverse locality in a dilute regime. 

Next question is about evolution equation for hardron wave function. In spirit of a 
QCD calculations at first approximation we can write result of infinitesimal boost as 

I _ / 1 + / G{z - x)a) {z)N {x)dxdz+ 
' ^ " + ] G*{z - x)G{z - y)N{x)N{y)dxdydz 

where universal function G{x) is a patron emission amplitude. It is directly related to the 
pomeron intercept and slope. This function must be extracted from more fundamental 
theory. Equation (j55p is the direct analog of the JIMWLK equation in the QCD. The 
pomeron here is just asymptotic of the scattering amplitude. 

In pomeron calculations widely used eikonal model with gauss vertexes 

iV(A:i...fc„) = <7"e-^'S'^ (56) 

It can be shown that in the parton model this equivalent to gauss distributions too 

/«(x)~e-'^^' (57) 

Unclear question here is about partonic interpretation of reggeon signature factor which 
is not depend on Y and on hardron properties. 
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